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We study the properties of a trapped interacting three component Fermi gas. We assume that one 
of the components can have a different mass from the other two. We calculate the different phases 
of the three component mixture and find a rich variety of different phases corresponding to different 
pairing channels, and simple ways of tuning the system from one phase to another. In particular, we 
predict co-existence of several different superfiuids in the trap, forming a shell structure, and phase 
transitions from this mixture of superfiuids to a single superfluid when the system parameters or 
temperature is varied. Such shell structures realize superfiuids with a non-trivial spatial topology 
and leave clear observable signatures in the density profile of the gas. 
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The observation of fermion pairing 0, Q, H, Q, Q 
and quantized vortices @ in ultracold Fermi gases have 
opened up new opportunities for exploring many-body 
quantum physics. For example, these systems can be 
used to study also strongly interacting fcrmionic su- 
perfiuids with mismatched Fermi surfaces Q- Experi- 
ments with such trapped polarized Fermi gases have re- 
vealed exciting and non-trivial phase separation proper- 
ties [1, 0, [T3| • These system have some intriguing simi- 
larities with electronic systems where external magnetic 
field imposes unequal mixture of electron spins as well 
as with quark matter systems where color superconduc- 
tivity may appear In color superconductivity mis- 



matched Fermi surfaces are due to unequal quark masses 
and such phenomena might be realized in the core of 
the compact star or in low-energy heavy- ion collisions. 
There is no fundamental problem in trapping and cool- 
ing also atoms with different masses and indeed both 
boson-boson as well as boson-fermion mixtures composed 
of different _atoms have been already experimentally re- 
alized 
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13l | and several groups are currently working 
towards the creation of different fermion-fermion mix- 
tures. The prospect of degenerate fermion mixtures with 
unequal masses has naturally motivated also recent the- 
oretical contributions 14, 15, 0, 17]. In this article 
we present general phase diagrams for a trapped three- 
component mixture of interacting fermions and predict, 
for instance, co-existence of several superfiuids and su- 
perfluid shell-structures. 

Some studies related to three-component Fermi gases 
have been done [H, 0| , but these studies were re- 
stricted either to an unlikely symmetric setup [Hj], to 
equal masses or additional simplifying assumptions 
were made concernin g th e strength of the various cou- 
plings in the system |l9l . [20j . Future experiments will 
certainly be done in an external trap and when the num- 
bers of particles are different, the trapped physics can be 
quite different from the physics in a uniform system. For 
this reason we will not here only allow different masses 



and varying coupling strengths in different pairing chan- 
nels, but will also solve the realized superfluid states in a 
trap. When interactions in all possible pairing channels 
are included, the number of trapped phases can become 
very high. For example, part of the trap might be oc- 
cupied by one type of superfluid, while some other part 
of the trap is occupied by a different superfluid or by a 
normal state. Even though the number of possible config- 
urations at first appears daunting, using the local density 
approximation, we find a simple picture which explains 
most of the trapped phases. As it turns out, one can ex- 
perimentally tune the system from one phase to another, 
by changing either the strength of interactions, trapping 
frequency, atom number, or temperature. 

We assume here that the third component has a differ- 
ent mass from the other two components. Consequently, 
as a by product, we also solve the trapped physics of the 
unequal mass two-component Fermi mixtures. Further- 
more, we briefly discuss the strength of induced interac- 
tion as the mass ratio is varied and find that the role 
of induced interactions becomes more important as the 
mass difference increases. 

We employ generalization of the BCS mean-field the- 
ory which is expected to be quantitatively fairly reliable 
in the weak coupling limit or with stronger interactions at 
zero temperature. At non-zero temperatures mean-field 
theory gives qualitatively reasonable results, but many 
body effects can alter the critical temperature substan- 
tially. Our Hamiltonian is given by 

H = /dr ^=1,2,3 (S + V*(r) - Av) &(r)Vv(r) + 
Ai 2 $(r)$(r) + A* 12 Mr)Mr) + Ai 3 $(r)$(r) 



A 23 ^(r)V4(r) 



|A 12 | 2 / 5 i 2 



|Ai 3 |7 ff i3 



+ A23^ 2 (r)4(r) 
|A 23 |7ff 3 2, (1) 



where V'i( r ) an d Vv(r) are fermionic creation and an- 
nihilation operators of the cr-component and interaction 
strengths are given by g aa ' = lixh 1 a aa i j '/x r , where a aa i 
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is the scattering length and fx r — (l/m a + l/m^) -1 is 
the reduced mass of the scattering atoms. For concrete- 
ness, most of our calculations are applied to a situations 
in which two of the components are different hypcrfinc 
states of 6 Li atoms and one of the components consists 
40 K atoms, although we also draw some general conclu- 
sions with arbitrary mass ratios. This implies mi = m-i- 
A ffff /(r) = g aa i (i/v(r)i/v(r)) are the order parameters 
which we find by minimizing the grand canonical poten- 
tial Q = -fc B TlnTr[e - ^]. Here the k B is Boltzmann 
constant, T is the temperature, and (3 — l/fc^T. Un- 
less indicated otherwise, we use the unit of length 1/kp, 
where kp is the Fermi momentum of the first component 
ideal gas in the center of the trap. Due to the use of 
contact interaction some integrals in momentum space 
diverge and are regularized, as usual, by subtracting the 
divergent contribution away. In this article we will focus 
on the most interesting case when all coupling strengths 
are negative. 

Here we take the trapping potentials to be harmonic 
V a (r) = m cr w^r 2 /2. The potentials are identical for the 
same mass atoms, which is accurate for the lowest two 
hyperfine states of 6 Li, for example, but the heavy atoms 
are assumed to experience a different trapping potential. 
This is in general true and having identical trapping po- 
tentials would require extra effort. In order for the details 
of our calculations to be valid, the heavy fermion should 
nevertheless experience a trap with a same aspect ratio as 
the light fermions. Most of our results will nevertheless 
remain qualitatively valid even with more complicated 
trapping scenarios. 

In order to compute different configurations we 

employ a local density approximation (LDA). In LDA 
we can, at each point in space, use a local chemical po- 
tential \i a (r) = fi a — V a (r) and find the state which min- 
imizes the grand potential. The chemical potentials in 
the center of the trap must then be determined in such a 
way that atom numbers have desired values. LDA could 
also be used in a more refined way by minimizing energy 
over the whole system [2l[ while including the possibil- 
ity of more exotic superfluid phases such as the breach- 
pairing/Sarma state. 

Trapped phases. — The minimization of the grand po- 
tential reveals that in a uniform system one has four 
different phases, three different superfluid phases corre- 
sponding to different pairing channels with the remaining 
component in a normal state, and a normal Fermi gas. 
We find some general rules in the weak coupling regime 
which predict the phase which is energetically favorable 
at zero temperature. First, pairing will occur only for 
closely matched Fermi surfaces. How closely the surfaces 
must be matched depends on the coupling strengths and 
stronger coupling implies greater tolerance for a Fermi 
surface mismatch. Second, if several pairing channels 
have same coupling strengths the pairing will occur in 
the channel with involving the heavy fermion. This ef- 



fect is due to the density of states increasing with mass 
which translates into lower energy. However, it should 
be kept in mind that this density of states effect can be 
masked by increasing the coupling in the channel involv- 
ing the light fermion. 

In the trap, the situation easily becomes more rich. 
The variables one can easily vary independently are the 
particle numbers of each component N a , ratio of the trap- 
ping frequencies for 40 K and 6 Li atoms, T, and using a 
(wide) Feshbach resonance, at least one of the coupling 
strengths. As we demonstrate in Fig. Q] at T = 0, one 
can easily find examples of systems with all four differ- 
ent phases co-existing in the trap at the same time. This 
is made possible by unequal atom numbers and different 
trapping frequency of the 40 K-component. 

The Fermi surfaces are matched when local chemical 
potentials satisfy /z<r(r) = m a jm a i\i a i{v). When the 
numbers of atoms in the equal mass 12-channel are differ- 
ent also their chemical potentials in the center of the trap 
are different and their Fermi surfaces are not sufficiently 
matched throughout the system. For a two component 
system this would result in a superfluid core surrounded 
by a normal gas @, 0, QjJ 0| , but in an interacting three 
component system superfluidity in the two unequal mass 
pairing channels can appear. Since the number of atoms 
and trapping frequencies are different, Fermi surfaces can 
be matched for different pairing channels at different lo- 
cations in the trapped gas. The requirement of matched 
Fermi surfaces implies that it is most natural to define a 
dimensionless effective trap frequency for the third com- 
ponent through r w = ?7iK^K/"^Li^Li- When r u = 1 the 
Fermi surfaces which are matched at the origin remain 
matched throughout the trap. For concreteness we as- 
sume a harmonic trap for Lithium atoms with a geomet- 
ric average of trap frequencies ui\a = 670Hz 

Since it is experimentally easier to observe integrated 
densities as opposed to a three-dimensional density dis- 
tributions, in part (b) of Fig. Q] we show the doubly in- 
tegrated density differences An aa i(z) — J dxdyn a (r) — 
n a i (r) . We find that superfluid regions appear as fiat ar- 
eas in the integrated density-differences. If the number of 
atoms participating in the superfluidity in some channel 
is high enough, it can leave behind a clearly observable 
experimental signature in the density distributions, not 
unlike what happens in a polarized two-component sys- 
tems num. 

Let us now discuss the dependence of the system on 
experimentally easily tunable parameters. In Fig. [2] (a) 
we show the phase diagram as a function of the effective 
frequency ratio r w and position in a trap. For equal effec- 
tive trapping the system is a pure 23-superfluid, due to 
the higher density of states for the unequal mass channel. 
The area occupied by the 23-superfluid becomes smaller 
with increasing and contracts to a shell surrounded 
by 12-superfmids. This is caused by the reduction in the 
area of sufficiently well matched Fermi surfaces with in- 
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creasing r w . 

Fig. [2] (b) shows the phase diagram as a function of 
the particle number ratio N3/N2. Because we choose 
r u = 1, /v(r) profiles have the same shape for all three 
components and one can get matched Fermi surfaces and 
superfluidity in the 23-channel only if the density of the 
third component is high enough in the center of the trap. 
For this reason there is no superfluid state A23 for small 
number of N3 atoms. Similar conclusion applies if the 
number of atoms in the third component is excessively 
high. 

In Fig. [3] (a) we show the phase diagram in the trap as 
a function — kpa\2- We find that when one approaches 
the Feshbach resonance and — kp a%2 becomes large there 
is phase transition from the superfluid state A23 to the 
superfluid state A12. In this case there are two com- 
peting two component systems. As the interaction in 
the same mass channel increases, pairing in that chan- 
nel eventually becomes energetically favorable even if the 
average mass in that channel is lower. In Fig. [3] (b) we 
demonstrate the phase diagram in the trap as a func- 
tion — kp ai2, when the trapping potentials are different. 
When kpayi S [—0.35,-0.2] we find, an unusual con- 
figuration, a normal gas in the center of the trap and 
a shell of 23-superfluid state surrounding it. This can 
be understood by the fact that pairing with the third 
component modifies the density distribution of the com- 
ponents so much that the Fermi surfaces of the first and 
second components do not match in the center of the trap 
anymore. A BCS solution is possible, at zero tempera- 
ture, only when — (12)/^ < A12 and since the gap 
A12 would be small for small I&FCI12I, modifications of 
density distributions do not have to be dramatic to make 
the BCS solution disappear in the center of the trap. 

The mean-field critical temperature for superfluid- 
ity between atoms with unequal masses in the BCS 
regime can be calculated and is given by ksT c — 
[S>ly/r^/{e 2 TT)]eFexp[TT/(2\kpa\)], where r rn = mu/rriK- 
Therefore, for the same value of the scattering length 
the T c for a system with unequal masses is changed by 
a factor of yjr m from the equal mass case. On the other 
hand at T — pairing in the unequal mass channel is 
energetically favored due to higher density of states of 
heavy fermions. Therefore, there is an intermediate (2nd 
order) transition from the unequal mass channel to the 
same mass channel. This is demonstrated in Fig. |4l 
where we show the phase diagram as of function of T. 
One clearly sees that as the temperature rises there is 
first a transition between the pairing channels, which is 
then followed by a transition into a normal state. In a 
trap, different pairing channels can again coexist in some 
intermediate temperature window. 

Many-body effects due to scattering in a medium mod- 
ify interactions between fermions and can substantially 
lower the critical temperature. This correction due to in- 
duced interactions was calculated by Gorkov and Melik- 




FIG. 1: (a) Gaps as a function of position (Rtf is the ideal 
gas Thomas- Fermi radius of the second component). Solid, 
dashed, and dashed-dotted lines describe A23, A13, and A12 
respectively, (b) Doubly integrated density differences. Solid, 
dashed, and dashed-dotted lines describe Ari23, An.13, and 
Ani2 respectively. We used the parameters iVi = 6-10 4 , N2 = 
5 ■ 10 4 , and N 3 = 2 • 10 4 , k F a 12 = -f.04, (f + l/r m )k F a 13 = 
— 1.03, (1 + l/r m )feFa23 = —1.15, where r m = m,Li/m,K and 
W3/W1 = 0.4. 
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FIG. 2: (a) Zero temperature phase diagram as a function r„ 
and t/Rtf with particle numbers Ni^2,3 = 5.5 ■ 10 4 . (b) Zero 
temperature phase diagram as a function particle number ra- 
tio N3/N2 and t/Rtf, with N\ = N 2 = 5.5 ■ 10 4 and = 1. 
In both Figures we used kFdi2 = —1, (1 + l/r m )kFa,23 = — 1, 
and |a23 1 > |cii3|.N indicates a normal gas or a vacuum at 
sufficiently large distances. 



Barkhudarov [23j and we have generalized this calcu- 
lation to deal with interactions between unequal mass 
fermions. We consider interacting atoms on the matched 
Fermi surfaces with energies equal to the chemical poten- 
tials since in the BCS regime physics close to the Fermi 
surface dominates. For the equal mass case the correc- 
tion in the critical temperature is the usual ~ 2.2, but 
it increases as the second component becomes heavier. 
With r m = 0.15 the correction turns out to be w 6.3. 
Due to this mass dependence of the Gorkov correction, 
we expect that finite temperature mean-field calculations 
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(a) (b) 




FIG. 3: Phase diagrams in the trap as functions —k F a\ 2 and 
t/Rtf for two different relative trap frequencies, in (a) r u ~ 
1.0 and (b) = 2.67. We used parameters 2Vi,2,3 = 5.5 • 10 4 , 
(1 + l/r m )k F a 2 3 = 1, and |<x 2 3 1 > |ai3 1- 



that in our case the surface contribution compared with 
the superfluid energy in the shell structures is at worst at 
the 2% level and usually much less than this. Therefore, 
they are unlikely to have a large effect on our phase di- 
agrams. Furthermore, the surface energy compared with 
the superfluid bulk energy scales in such a way that it can 
always be made small by either increasing atom numbers 
or by reducing the trap frequency. It would be inter- 
esting to understand how co-existing superfluids of non- 
trivial spatial topology influence the collective and the 
rotational properties of these systems. 

This work was supported by Academy of Finland 
(project numbers 207083, 106299, 205470) and con- 
ducted as part of a EURYI scheme award. See 
www.esf.org/euryi. 
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FIG. 4: Phase diagram as a function of the temperature T 
and position. We used parameters Ni,2,3 = 5.5 • 10 4 , — 
1, k F a 12 = -1, (1 + l/r m )k F a 2 3 = -1, and |a 2 3| > |ois|. 
With our parameters the Fermi temperature T F ~ 600nK 
and Rtf ~ 50 /im. 

overestimate the critical temperature in the channel with 
unequal masses relative to the equal mass channel. 

Conclusions. — We have solved the mean-field theory of 
a trapped three-component Fermi gas We found different 
superfluids coexisting in a same trapped system which 
can leave clear signals into the easily measurable density 
distributions. Such signatures of the spatial distribution 
of the different superfluids could be complemented by in- 
formation about the different pairing gaps obtainable by 
rf-spectroscopy [B|, [24[ ■ With the our example parame- 
ters the coherence length in the center of the cloud turns 
out to be 0.11 • Rtf- While this is a fairly large fraction 
we nevertheless choose to use the local density approxi- 
mation. This is unlikely to affect the general features of 
our results, which hinge on pairing with matched Fermi 
surfaces. By increasing the atom numbers the local den- 
sity approximation becomes ever more accurate In this 
article we have ignored the surface energy contributions 
between different superfluids. However, by applying the 
results of De Silva and Mueller [25J we have estimated 
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